Nonlinear screening of charged 
macromolecules 



By Gabriel Tellez 
Departamento de Fisica, Universidad de los Andes, A. A. 4976, Bogota, Colombia 

We present several aspects of the screening of charged macromolecules in an elec- 
trolyte. After a review of the basic mean field approach, based on the linear Debye- 
Hiickel theory, we consider the case of highly charged macromolecules, where the 
linear approximation breaks down and the system is described by full nonlinear 
Poisson-Boltzmann equation. Some analytical results for this nonlinear equation 
give some interesting insight on physical phenomena like the charge renormaliza- 
tion and the Manning counterion condensation. 
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1. Introduction and the linear Poisson— Boltzmann equation 

A colloidal suspension is a system composed of two substances, one dispersed into 
the other. The dispersed phase is composed of macromolecules, with size of the 
order 10~* m - 10~^ m, while the dispersion medium, or continuous medium, is 
composed of small micromolecules, and/or microions, with size of the order of the 
nanometer. Of particular interest are the charge-stabilized colloids, where the dis- 
persed macromolecules have ionizable sites and, when immersed into the dispersion 
medium, they acquire a surface electric charge which ensures repulsion between 
them and thus allows the colloid to stabilize and prevents aggregation. 

Since the length and time scales of the microcomponents of the dispersion 
medium are much smaller than the ones of the dispersed phase, it is convenient 
to average over the degrees of freedom of the dispersion medium and treat the 
system as a one-component system composed by the macromolecules which inter- 
act via an effective potential. To understand the physical and thermodynamical 
properties of these systems, it is important to determine this effective interaction 
between the macromolecules, which results not only from the direct interaction be- 
tween the macromolecules, but also the interaction mediated by the microions of 
the dispersion medium. 

The basic theory to find this effective interaction for charge-stabilized colloid was 
developed independently by [Derjaguin and Landau 1941] , and Overbeek and Verwey 1948 



and it is known as the DLVO theory. It is based on the work of [Debye and Hiickel 1923| 
Let us consider a spherical charged macromolecule, with radius a, and charge Ze 
(e is the elementary charge) immersed in an electrolyte with positive microions of 
charge z+e and average density n_|_, and negative microions of charge — z_e and aver- 
age density n_ . Without loss of generality we can suppose Z > 0. [Debye and Hiickel 1923| 
idea is to treat the microions in a mean field approximation: at temperature T, the 
local density of microions of charge q± = ±z±e at a distance r from the macro- 



Article submitted to Royal Society 



Paper 



2 



Gabriel Tellez 



molecule can be approximated by 

n±(r) =n±e-'3«±*W, (1.1) 

where ^{r) is the electrostatic potential, and (3 = l/{kBT), with ks the Boltzmann 
constant. Replacing this into Poisson equation of electrostatics yield the Poisson- 
Boltzmann equation 

A* = -— (z+n+e-^^'+^ ~ z^n^e"^'-'^) , (1.2) 

where e is the dielectric constant of the dispersion medium. 

It is convenient to introduce the following notations: the reduced potential y = 
/3e^, the Bjerrum length Ib = /3e^/e, the Debye length — (47rZB(z^n+ + 
z^n_))^^/^. With these notations, Poisson-Boltzmann equation reads 

Ay = \e'-y - e-'-+y] . (1.3) 

z+ + Z- 

If the electrostatic coupling between the macromolecule and the microions is small, 
y{r) <^ 1, for any distance r, the nonlinear Poisson-Boltzmann equation p. 31) can 
be linearized to obtain 

Ay^K'y. (1.4) 

For an impenetrable spherical macromolecule with uniform surface charge (total 
charge Ze and radius a), the solution of this equation is the DLVO potential 

y{r) = ZlB— . (1.5) 

1 + Ka r 

From this equation, one can sec that 1/k is the screening length. 

It is also interesting to consider the case of cylindrical macromolecules, for in- 
stance stiff polyelectrolytes, ADN, etc. As a first approximation, for an infinitely 
long cylinder with radius a and linear charge density e/£, uniformly spread over its 
surface, the solution of equation (|1.4p gives the electrostatic potential at a radial 
distance r from the cylinder 

y{r) = ^r^^Ko{f), (1.6) 

where Kq and Ki are the modified Bessel functions of order and 1 . We have defined 
the reduced linear charge density of the cylinder ^ = and it is convenient to 

measure the distances in Debye length units: f = nr and a = na. The boundary 
conditions that complement the differential equation (|1.4p to yield the solution (|1.6|) 
are 

lim = -2C, and limV?;(r) = 0. (1.7) 

It should be noticed that at large distances from the cylinder, compared to the 
Debye length, the potential exhibits again an exponential decay, as for the case of 
spherical macromolecules. 



yM-.r^W;5^e— , r » . (1.8) 
aKi[a) V ZKr 
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In the following sections, we turn our attention to the case of highly charged macro- 
molecules, where the linear approximation breaks down and the full nonlinear equa- 
tion (|1.3p should be used. Several nonlinear effects appear, which we will review. 
In section [5J we study the nonlinear phenomenon known as the charge renormal- 
ization [Alexander et al. 1984[ ITrizac et al. 2002] . which is generic, both for the 
spherical and cylindrical geometries. In section [31 we focus our attention on the 
cylindrical case, where an analytical solution of the nonlinear Poisson-Boltzmann 
equation is available. From the analysis of this analytic solution we discuss the 
phenomenon of countcrion condensation. 



2. Charge renormalization 

For a highly charged macromolecule, Poisson-Boltzmann equation ([1.3p cannot be 
linearized near the macromolecule surface. However, due to the screening effect, 
the potential will decay and become small, \y{r)\ <^ 1, at large distances from the 
macromolecule surface, r — a ^ . In that far region, the linear version ([1.4p 
of Poisson-Boltzmann equation holds. Then, for a spherical macromolecule, the 
potential will behave as 

y{r) ^ a"-^ (2.1) 

at large distances from the macromolecule. To find the constant of integration 
one needs to enforce the boundary condition at the surface of the macromolecule 
that the normal component of the electric field is proportional to the surface charge 
density. However, the form (|2.ip of the potential is not valid in that close region. 
One needs to find also the form of the potential close to the macromolecule sur- 
face, and connect it to the large- distance behavior (|2.ip to find explicitly the in- 
tegration constant A. In analogy to the linear solution (|1.5p . one can write the 
A = Zj-cuIbg'^"' /{I + Ka), defining a renormalized charge Zrcn- The large-distance 
behavior of the potential then takes a DLVO familiar form 

y{r) ^ Z,,Jb — , (2.2) 

1 + Ka r 

but replacing the bare charge Z of the macromolecule by the renormalized one .^ren- 
in the cylindrical geometry, the renormalized charge concept also applies. In that 
case, the large-distance behavior of the potential is 

yW-Tf^ifo(r), (2.3) 

with a renormalized linear charge density ^ron- The determination of the renormal- 
ized charge requires knowledge of the short-distance behavior of the solution of 
the nonlinear Poisson-Boltzmann equation ([1.3p . This can be done numerically, as 
in the original work of [Alexander et al. 1984j . In experimental situations, the un- 
known renormalized charge is often taken as an adjustable fitting parameter. There 
are also analytical approaches to find the renormalized charge [Trizac et al. 20021 
IShkel et al. 20001 ITellez and Trizac 2004j . mostly based on approximations using 
the solution to the nonlinear Poisson-Boltzmann equation in the planar case. 

Let us illustrate the concept of charge renormalization in the planar geome- 
try where an explicit solution for the nonlinear Poisson-Boltzmann equation is 
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known |Gouy 1910[ Chapman 1913| . The system is an infinite charged plane, with 
charge density cr > 0, immersed in an electrolyte, which, for simplicity, we consider 
symmetric = z_ = 1. Let Ox be the axis perpendicular to the plane, which we 
suppose located at x = 0, the electrolyte occupies the region a; > 0. The nonlinear 
Poisson-Boltzmann equation in this situation reads 

V ' = sinh y (x) . (2.4) 
ax 

It can be integrated once by multiplying by dy/dx, 

M^ = _2.sinh^. (2.5) 
dx 2 

where the boundary condition dy/dx —^ when x ^ cxi has been used. This last 
equation is separable and can be integrated, finally obtaining 

2/W = 21n ^_^^_.. . (2.6) 

Where A is a constant of integration, which is found using the boundary condition 
at the surface of the charged plane 

^(0) = -47rWe- (2.7) 
dx 

Notice that at large distances from the plane, kx ^ 1, the potential behaves as 

y{x) ~ 4Ae-'^^ . (2.8) 

This is the expected behavior for the linear version of Poisson-Boltzmann equa- 
tion in this geometry y{{^{x) — K^yiin(x) — 0. The linear solution is yiin(2;) — 
AnlBcre^'^^ / [en). Comparing to the large-distance behavior of the nonlinear so- 
lution (j2.8p . one can define the renormalized surface charge density (Trenj by writing 
the constant of integration A as 

A = TrlB<Jrcn/{eK) . (2.9) 

Then, the large-distance behavior of the potential is 

y{x) e (2.10) 

To find explicitly the renormalized charge, one needs to apply the boundary condi- 
tion (|2.7p at the surface of the charged plane. For this, one needs the short-distance 
behavior of the potential. From the explicit solution (|2.6p . we find 

1 + A AAx 
1-A ^ 1-^2 



y(x) = 21n^ -^ + 0(x2), (2.11) 



when X ^ 0. Using this, we apply the boundary condition (|2.7p to find 

A Train 
1 J?^ -■ 2.12 
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Figure 1. Renormalized charge as a function of the bare charge in the planar case, for a 

symmetric electrolyte. 

Solving, and using (|2.9p . we find the renormalized surface charge 

a 

where we have defined reduced charge densities (7 — -kIbo j (en) and (Tj-cn — '^Is'^rcn/ 
Figure [1] shows a plot of the renormalized charge as a function of the bare charge. 
Notice the saturation effect: when a cxi, the renormalized charge approaches 
a finite value (Jren — *■ 1- This saturation effect can also appear in other theories 
obtained by modification of Poisson-Boltzmann equation using a density functional 
formulation [Tellez and Trizac 2003] . In that saturation regime, the large-distance 
behavior potential becomes independent of the bare charge of the plane 

ysat(x)~2;oe~"" = 4e-^^. (2.14) 

The value yo = 4 plays an important role. It can be seen as an effective surface po- 
tential for the plane, if one wants to match the linear solution of Poisson-Boltzmann 
equation with the nonlinear one in the close vicinity of the highly charged plane. 
This is also the starting point to find the renormalized charge at saturation for 
highly charged macromolecules of arbitrary shape. For these macromolecules, when 
they are highly charged, the linear Poisson-Boltzmann equation can be solved with 
an effective boundary condition of constant surface potential yo, to find the behav- 
ior of the potential at large distances and the corresponding renormalized charge, 
as explained in [Bocquet et al. 2002| . For instance, for a spherical macromolecule, 
the solution of the linear Poisson-Boltzmann equation with the effective constant 
potential boundary condition at the surface of the macromolecule y{a) = j/o = 4, is 

2/sat, sphoro('') = yoc'^^O . (2.15) 



The nonlinear solution, at large distances has the behavior given by equation 
Comparing both equations (|2.2p and (|2.15p . we find an approximate value for the 
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renormalized charge in the saturation regime 



Zrcn = -^^yoina + 1) = -^{Ana + 4) . 



(2.16) 



This approximation is based on the planar solution of the nonlinear Poisson- 
Boltzmann equation, and therefore it is accurate for large macromolecules with 
Kfl 3> 1, and only at the first order in Ka. One can improve this estimate by 
developing a planar expansion of the solution of the spherical geometry, as done 
in [Shkel et al. 20001 ITrizac et al. 20031 ITellez and Trizac 2004] . Up to terms of or- 
der 0(l/(«;a)), the renormalized charge at saturation for spheres in a 1:1 electrolyte 
is [Trizac et al. 2003] 



h 



-(4Ka + 6), 



(2.17) 



The starting point to obtain estimates of the renormalized charge at saturation 
is the value yo of the effective surface potential at saturation in the planar geometry. 
This depends only on the constitution of the electrolyte. We propose now a simple 
formula that gives yo in the generic case of a multicomponent electrolyte, composed 
of several species of ions with charges {qae} and densities {na}. The nonlinear 
Poisson-Boltzmann equation in the planar geometry reads now 



Ha 



„e^9°«(^) = 0. 



(2.18) 



Multiplying by y'{x) this equation, it can be integrated once to find 

(y'(x))2 =871^5 ^ (e-«°«(^)-l), (2.19) 



where the boundary condition y'{x) when x — > oo has been used. Introducing 
the inverse Debye length k = {AttIb X)q 'Za'^a)"^^^, one obtains the formal solution 

KX = / . (2.20) 



Suppose the charged plane is located at x = and positively charged, and we are 
in the saturation regime, therefore j/(0) +00. At large distances from the plane, 
KX ^ 1, the potential behaves as 7/(a;) ~ ?/oe~'^'^, thus = \n.yQ—\n.y{x)+o{\n.y{x)) 
as y{x) — > 0. Replacing in (|2.20p we find 



Inyo 



lim 



= Pf. 



du 



du 



Iny 



(2.21) 



Thus, the value of j/o is expressed as a Hadamard finite part (Pf.) of the inte- 
gral (|2.2ip . Alternatively, it can be computed from 



Inyo = 



1 



— E."a(e-''""-l) 



du' 



du 



(2.22) 
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Table 1. Value of the saturation potential for several electrolytes. 



2+ : Z- 


1:5 


1:4 


1:3 


1:2 


1:1 


2:1 


3:1 


4:1 


5:1 


yo 


0.56080 


0.71742 


0.99388 


1.6077 


4 


6 


8.7070 


12.314 


17.337 






z+ : Z- 


2:3 


3:2 


2:5 


5:2 




3:5 


5:3 


1:10 


10:1 


yo 


1.1542 


2.4611 


0.61471 


3.6270 




0.67242 


1.7544 


0.26761 


70.337 



In the case of a two-component electrolyte qi = z+ and 52 = —z-, the saturation 
value j/o can be expressed as a function of the ratio r = z+/ z-, 

1 

I/O = — exp 



In the cases of electrolytes with z+ : z_ equal to 1:1, 1:2, and 2:1, the integral can 
be computed exactly to find the known values [Tellez and Trizac 2004) . y}f-^ = 4, 
— 6(2 — \/3) ~ 1.6077, and j/q'^ = 6. Table [T] gives the value of the saturation 
potential for other electrolytes. 




(2.23) 



3. Cylindrical macromolecules 

(a) Exact solutions for Poisson-Boltzmann equation and the connection problem 

In this section, we focus our attention on the study of the screening of a 
thin cylindrical macromolecule, with radius a ^ k^^. Besides the charge renor- 
malization effect, another interesting phenomenon that occurs in this geometry 
is the counterion condensation. This was first realized by Onsager and studied 
by [Manning 1969] and [Oosawa 1971] . To understand this phenomenon, consider 
the Boltzmann factor between the macromolecule and an ion of opposite charge 
(counterion): exp(— 2z_^lnr) = 7--22_c if^ diverges when r ^ 0, and furthermore 
it is not integrable near r — > if ^ > 1/z-. This means that, for an infinitely thin 
macromolecule, with radius a — 0, the thermodynamics are not properly defined 
unless ^ < l/z„. In real situations a 7^ 0. For ko ^ 1, when ^ > l/z-, the density 
of counterions will be very large near the surface of the macromolecule. These coun- 
terions are bounded to the macromolecule: besides the diffuse screening cloud of 
ions around the macromolecule, there is also a thin layer of condensed counterions 
very near to the surface of the cylinder. 

This counterion condensation effect can be studied quantitatively in the mean 
field approximation, since an analytical solution of the nonlinear Poisson-Boltzmann 
equation in the cylindrical geometry is available. For a 1:1 electrolyte, this so- 
lution was found by [McCoy, Tracy and Wu 1977 , in a different context, in re- 



lation to the correlation functions of the two-dimensional Ising model. Later on, 
[Widom 1997| developed the solution for the asymmetric cases 2:1 and 1:2, and 
[Tracy and WidonTlQ QS studied the short-distance asymptotics of the solution and 
solved the problem of connecting the large-distance and the short- distance behav- 
iors of the solution. The consequences of this mathematical work to the screen- 
ing of cylindrical macromolecules were reported by [McCaskill and Fackerell 1988j , 
[Tracy and Widom 1997[ and [Trizac and Tellez 2006| . A more extensive study is 
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presented in [Tellez and Trizac 2006] . We summarize here some of the main findings 
of that work. 

The short-distance behavior of the potential can be obtained by a physical argu- 
ment. Very close to the charged cylindrical macromolecule, one would expect the po- 
tential to be the bare Coulomb potential — In r -I- constant, with A some constant 
related to the charge of the macromolecule. By replacing this ansatz into Poisson- 
Boltzmann equation (|1.3p one can compute systematically the following terms of the 
short-distance expansion, to find [Tracy and Widom 1998| ITellez and Trizac 20061 
ITrizac and Tellez 2007| 



y{r) = -2Alnf + 21nB - 2 In 



2 



1 - 



16(1-^)2 



0( 



s2+2An 



(3.1) 



We specialize here in the 1:1 electrolyte. For the general case of a z+:z- electrolyte, 
see [Trizac and Tellez 2007[ . A and B are some constants of integration. The con- 
stant A can be related to the charge density of the electrolyte by writing the first 
boundary condition (|1.7p 



(2 - 2A){Ka) 



2-2A 



B^ 



16(l-A)2-B2(Ka)2-2A 



(3.2) 



From the previous section, we already know the large-distance behavior of the 
potential, it is the screened potential (|2.3p 



y{r) ^ AXKaif) , 



(3.3) 



where A is some constant related to the renormalized charge by 



6cn = aifi(a)A/2. 



(3.4) 



By using the explicit analytical solution of the nonlinear Poisson-Boltzmann equa- 
tion from [McCoy, Tracy and Wu 1977] and ^Widom 1997| . [Tracy and Widom 1998] 
where able to solve the connection problem of relating the constants of integration 
from the short-distance behavior A and B, to the one of the large-distance behavior 
A. To satisfy the boundary condition y'{r) when r oo, the constants A and 
B need to satisfy 

B^2^^ ^\/\i , (3.5) 



' l+A -\ 
- 2 ) 

V 2 ; 



where F is the gamma function, and, A and A need to satisfy 



1 . f7:A 
A = — sm — — 
TT \ 2 



(3.6) 



The first physical consequence of these relations, is that we can obtain an ana- 
lytical expression for the renormalized charge, by combining equations (|3.2[) . (|3.4p 
and ()3.6p . In the simplest situation, when a = 0, this gives, A — ^ and 



, 2 . TT^ 

4ron = - sm — 
TT I 



(3.7) 
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ih) Counterion condensation 

The previous discussion, and in particular equation p.ip are only valid provided 
that A < 1. Indeed, if A — 1, equation (|3.1[) becomes singular: the last term becomes 
of the same order as the second, and besides, the constant B from equation ()3.5p 
becomes undefined. This is the mathematical signature of the counterion conden- 
sation phenomenon. Notice that for a = 0, the constant A is the linear charge of 
cylinder A = ^, and the value ^Manning = 1 is precisely the threshold for counterion 
condensation discussed earlier. 

For a cylinder with nonzero radius a 7^ 0, notice that, using equation (j3.2|) . the 
threshold A — ^Manning = 1 coircsponds for the linear charge to the threshold 

?c = 1 + ^-r^ , (3.8) 
In a + O 

with C = 7 — 3 In 2 ~ —1.502, with 7 the Euler constant. Note that there is a 
negative logarithmic correction in the radius a of the cylinder to the Manning value 
^Manning = 1 for the threshold for condensation: ^ 1- 

To extend the solution beyond the condensation threshold, [Tracy and Widom 1998| 
suggested to write A as a complex number A = 1 — i^/2. Replacing into equa- 
tion p.ip . the short-distance expansion of the potential now reads 

, , „, sin(— 2/ilnf — 2uC) , . 

yfr = -2 Inf- 2 In — i ^-^ . (3.9) 

The constant fi can be expressed in terms of the bare linear charge density ^ by using 
the first boundary condition (|1.7p . and connected to the large-distance expansion 
of the potential and the renormalized charge by means of equation (j3.6p . replacing 
A — 1 — i/i/2, for details see [Tellez and Trizac 2006] . 

Notice the first term of expansion (|3.9p : — 2 In f. It is the bare Coulomb potential 
of a charged line with linear charge density ^Manning = 1- This a characteristic of 
the counterion condensation phenomenon. At intermediate distances of the charged 
cylinder, one "sees" a cylinder with an effective charge ^Manning ~ 1, if the bare 
charge exceeds the threshold value ^c- The second term of equation (|3.9p can become 
very large in the close proximity of the cylinder. This second term represents the 
thin layer of condensed counterions located at the surface of the cylinder. Figure [H 
shows a plot of the potential close to the cylinder, and compares it to the bare term 
-2 Inf. 



4. Concluding remarks 

To summarize and conclude the review presented here, we would like to stress 
the differences between the different linear charges densities that we presented for 
cylindrical macromolecules: renormalized charge. Manning charge, threshold charge 
for counterion condensation. 

The renormalized charge characterizes the behavior of the potential far from the 
charged macromolecule. At those large distances, r , the potential exhibits 

an exponential decay. The prefactor of this exponential decay is proportional to the 
renormalized charge, as shown in equation (j2.3p . 
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Figure 2. Short-distance expansion of the potential (|3.9|) (thick line) compared to 
— 21n(f/(4/i)) (dashed hne). Parameter n = 0.2. 

At short distances, there are two possible behaviors, depending on the value of 
the bare linear charge density ^ compared to the threshold value given by equa- 
tion (|3.8p . If < Cc7 the potential behaves as given by equation (|3.ip . It is a bare 
Coulomb potential with a prefactor given in terms the bare charge of the macro- 
molecule. If ^ > £,c, the counterion condensation takes place. The short-distance 
behavior of the potential is given by equation (|3.9p . A thin layer of counterions 
is bound to the surface of the cylinder, which makes the potential very large in 
that region. Beyond this layer, the potential behaves as — 2^Manninglnr : the bare 
Coulomb potential for a charged cylinder but with an effective charge ^Manning = Ij 
the Manning charge for counterion condensation, see figure [H Notice that if the 
radius of the cylinder a ^ 0, the Manning value differs from the threshold value: 

^Manning ^ ^c- 

The counterion condensation phenomenon can be only noticed at close proximity 
of the charged cylinder, by the change of the short-distance behavior of the poten- 
tial from (|3.1|) to (|3.9p . At large distances, the potential is always given by (j2.3|) . 
regardless if the counterion condensation has taken place or not. When ^ = ^c, no 
singularity appears in equation (j2.3p . nor in the renormalized charge ^ron which 
characterizes only the large-distance behavior of the potential. Also, notice that 
iron 7^ ^Manning and ^ren 7^ Cc- This has causcd somc confusion in the past, since in 
the original work of [Manning 1969| , in the condensed phase, the diffuse screening 
cloud of the remaining uncondensed counterions around the charged cylinder was 
treated using the linear Poisson-Boltzmann equation and using an effective charge 
of the cylinder given by ^Manning = 1 to account for the counterion condensation. 
It turns out that this picture is not completely correct, since besides the counte- 
rion condensation, there are also additional nonlinear effects in the uncondensed 
cloud that are responsible for an additional charge renormalization. As a result 
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Cicn 7^ ^Manning- For example, in the limiting case a = 0, from equation p.7p . we 

obtain ^rcn = S/tT ~ 0.63662 < ^Manning = 1- 

The author acknowledges partial financial support horn Comite de Investigaciones y 
Posgrados, Facultad de Ciencias, Universidad de los Andes. 
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